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Q^ , Abstract. We introduce a discrete-time fractional calculus of variations. First 

l/^ ' and second order necessary optimality conditions are established. Examples 

C 3 ' illustrating the use of the new Euler-Lagrange and Legendre type conditions are 

given. They show that the solutions of the fractional problems coincide with 
the solutions of the corresponding non-fractional variational problems when 
the order of the discrete derivatives is an integer value. 



1. Introduction. The Fractional Calculus is currently a very important research 
field in several different areas: physics (including classical and quantum mechanics 
. , and thermodynamics), chemistry, biology, economics and control theory [20, 23, 

j^ ■ 24, 27, 28]. It has origin more than 300 years ago when L'Hopital asked Leib- 

niz what should be the meaning of a derivative of order 1/2. After that episode 
several more famous mathematicians contributed to the development of Fractional 
Calculus: Abel, Fourier, Liouville, Riemann, Riesz, just to mention a few names. 

In [22] Miller and Ross define a fractional sum of order v > via the solution 
of a linear difference equation. Namely, they present it as (see Section 2 for the 
notations used here) 



A-V(i) = ^ E(^ - ^(■^)Y''-''^fi-^)- (1) 
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This was done in analogy with the Riemann-Liouville fractional integral of order 



aD-'-fix) = ^ J\x sr-^fis)ds, 



which can be obtained via the solution of a linear differential equation [22, 23]. 
Some basic properties of the sum in (1) were obtained in [22]. More recently, F. Atici 
and P. Eloe [9, 10] defined the fractional difference of order a > 0, i.e., A"/(i) = 
A'"(A~(™~")/(i)) with m the least integer satisfying m > a, and developed some of 
its properties that allow to obtain solutions of certain fractional difference equations. 

Fractional differential calculus has been widely developed in the past few decades 
due mainly to its demonstrated applications in various fields of science and engineer- 
ing. The study of fractional problems of the Calculus of Variations and respective 
Euler-Lagrange equations is a fairly recent issue - see [1, 2, 5, 7, 8, 11, 13, 15, 16, 21] 
and references therein - and include only the continuous case. It is well known that 
discrete analogues of differential equations can be very useful in applications [18, 19]. 
Therefore, we consider pertinent to start here a fractional discrete-time theory of 
the calculus of variations. 

Our objective is two-fold. On one hand we proceed to develop the theory of 
fractional difference calculus^ namely, we introduce the concept of left and right 
fractional sum/difference (cf. Definitions 2.1 and 2.4 below) and prove some new 
results related to them. On the other hand, we believe that the present work 
will potentiate research not only in the fractional calculus of variations but also in 
solving fractional difference equations, specifically, fractional equations in which left 
and right fractional differences appear. 

Because the theory of fractional difference calculus is in its infancy [9, 10, 22], the 
paper is self contained. We begin, in Section 2, to give the definitions and results 
needed throughout. In Section 3 we present and prove the new results; in Section 4 
we give some examples. Finally, in Section 5 we mention the main conclusions of 
the paper, and some possible extensions and open questions. Computer code done 
in the Computer Algebra System Maxima is given in Appendix. 

2. Preliminaries. We begin by introducing some notation used throughout. Let 
a be an arbitrary real number and 6 = fc + a for a certain /c G N with fc > 2. We put 
T = {a,a-hl,...,6},T'= = {a, a-h 1, . . . , 6- 1} and T'^' = {a, ...,&- 2}. Denote by 
F the set of all real valued functions defined on T. Also, we will frequently write 
(j{t) = i + 1, p{t) = i-1 and fit) = f{cr{t)). The usual conventions J2t=l fW = 0' 
c G T, and J^^q f{i) = 1 remain valid here. 

As usual, the forward difference is defined by A/(t) = f^it) — f{t). If we have 
a function / of two variables, /(i, s), its partial (difference) derivatives are denoted 
by At and A^, respectively. For arbitrary x,j/ S M define (when it makes sense) 

^iy) = r(x + i) 

r(x + l-y)' 
where F is the gamma function. The following property of the gamma function, 

F(x + l)=a;F(x), (2) 

will be frequently used. 

As was mentioned in Section 1, equality (1) was introduced in [22] as the frac- 
tional sum of order z^ > 0. While reaching the proof of Theorem 3.2 we actually 
"find" the definition of left and right fractional sum: 



FRACTIONAL DIFFERENCE PROBLEMS OF THE CALCULUS OF VARIATIONS 3 

Definition 2.1. Let f € J-. The left fractional sum and the right fractional sum 
of order i/ > are defined, respectively, as 

,Ar/W = ^E(^-'^G^))^^"'^/W> (3) 

^ ■' s—a 

and 

*Ar/W = f^ E (^ - a{t)Y''-'^f{s). (4) 

Remark 1. The above sums (3) and (4) are defined for t G {a+u, a+i/+l, . . . , b+v} 
and t G {a — u,a — ly + 1, . . . ,b — v}, respectively, while f{t) is defined for t e 
{a,a+l, . . . ,b}. Throughout we will write (3) and (4), respectively, in the following 
way: 

aA^fit) = J- ^(f + ^ - a(.))(''-^)/(.), t e T, 

^ ■' s—a 

Remark 2. The left fractional sum defined in (3) coincides with the fractional sum 
defined in [22] (see also (1)). The analogy of (3) and (4) with the Riemann-Liouville 
left and right fractional integrals of order j/ > is clear: 



,D-V(x) = -i^ /'(a 



X - sy-\f{s)ds, 



,D-v(^) = f^ / (^ - xr-'f{s)ds 



r(^) y. 

It was proved in [22] that lim^^o aA^'^ fit) ~ /(*)• We do the same for the right 
fractional sum using a different method. Let j/ > be arbitrary. Then, 

tA^-'fit) = f^ E('^ + ^ - <^it)Y''-'^fi^) 

^/W + T^T E {s + i^-<j{t)r-'^f{s) 

Therefore, lim^_j.o t^^'^ f{t) ~ f{t). It is now natural to define 

,A°/(i) = ,A°/(f) = /(t), (5) 

which we do here, and to write 

,Ar/(i)=/W + ^;^^;^E(^ + ^-^(^))*'"''/W' ^e^' ^^0' (6) 
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0. 



.a{t) 



The next theorem was proved in [10]. 
Theorem 2.2 ([10]). Let / G -F and u > 0. Then, the equality 

aA-^Afit) = AiaA-^fit)) li±iiZ^)^/(a), t e T« 



r(^) 



holds. 



Remark 3. It is easy to include the case i^ = in Theorem 2.2. Indeed, in view of 
(2) and (5), we get 



aA^Afit) = AiaA^fit)) 



(t + ^_a)(-i)/(a), teT^ (7) 



r(^ + i) 

for aU 1/ > 0. 

Now, we prove the counterpart of Theorem 2.2 for the right fractional sum. 
Theorem 2.3. Let J '^ T and v >0. Then, the equality 



^:aAm 



t^pib) 

holds. 



r(^ + i) 



(6 + i/-f7(<))(''-iV(&) + A(tAr/(t)), teT^ (8) 



Proof. We only prove the case ly > a,s the case i^ = is trivial (see Remark 3). 
We start by fixing an arbitrary t e T**. Then, we have that, for all s £ T'^, 

A,((s + z.-a(i))(''-^V(s))) 

^{u- l){s + V - a{t))^''-'^ris) + {s + u- a(t))(''-i)A/(s), 

hence, 



6-1 



^\2i-' + u-a{t))^-^-'^Af{s) 

^ ■' S — t 



(S + l^ - (7(t))('^-l) 



r(^) 



/(«) 



fc-i 



=t ^ ) s=t 



(^+--(^»'^'-^/(6)-(-;)r\(t) 



r(^) 



r(^) 



1 ^"' 



n-) ^, 
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We now compute A^A^" f{t)): 



Ht^riii)) 



1 



r(^) 



r(^) 



s=a{t) 
b 

s=t 
b 



i=a(t) 



(!.-l)('^-l) 



1 



Y At{s + u-n{t)) 



(1.-1) 



/(«) 






fit) 



fit) 



^ ' s=a{t) ^ ' 

^ Y^v^ i)(.s + V a(t))(-2)r(.) - ^" ,;^ fit). 



rwet 



r(i.) 



Since i is arbitrary, the theorem is proved. D 

Definition 2.4. Let < a < 1 and set [i = 1 — a. Then, the left fractional 
difference and the right fractional difference of order a of a function / G J-" are 
defined, respectively, by 

aAtfit) = AiaAf^ fit)), t e T\ 

and 

tA^fit) = -AitA-^fit))), t e T^ 

3. Main Results. Our aim is to introduce the discrete-time fractional problem of 
the calculus of variations and to prove corresponding necessary optimality condi- 
tions. In order to obtain an analogue of the Euler-Lagrange equation (cf. Theo- 
rem 3.5) we first prove a fractional formula of summation by parts. The results of 
the paper give discrete analogues to the fractional Riemann-Liouville results avail- 
able in the literature: Theorem 3.2 is the discrete analog of fractional integration by 
parts [26, 27]; Theorem 3.5 is the discrete analog of the fractional Euler-Lagrange 
equation of Agrawal [1, Theorem 1]; the natural boundary conditions (22) and (23) 
arc the discrete fractional analogues of the transversality conditions in [3, 21]. How- 
ever, to the best of the authors knowledge, no counterpart to our Theorem 3.6 exists 
in the literature of continuous fractional variational problems. 

3.1. Fractional Summation by Parts. The next lemma is used in the proof of 
Theorem 3.2. 

Lemma 3.1. Let f and h be two functions defined on T'^ and g a function defined 
on T'^ X T*^ . Then, the equality 

b-l T-l 6-2 6-1 

5^/(r)^5(r,s)M.) = EMr) ^ 9is,r)fis) 

T=a s=a T=a s=a{T) 

holds. 
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Proof. Choose T = Z and F(t, s) = .f{T)g{T, s)h{s) in Theorem 10 of [6]. 



D 



The next resuh gives a fractional summation by parts formula. 



Theorem 3.2 (Fraetional summation by parts). Let f and g be real valued functions 
defined on T'' and T, respectively. Fix < a < 1 and put /i = 1 — a. Then, 



6-1 



b-2 



Y^ f{t)aA'^g{t) = f{b l)g{b) f{a)g{a) + ^ ,A^(,)/(t).g'^(t) 



t=a 



h-l 



t — a 



\l^~^ ) \ i^a t=a{a) 



Proof. From (7) we can write 



6-1 



6-1 



5^ /(t)aA?5W - E /(OA(aAr^ffW) 
t^a t—a 

b-1 

= E/w 



t — a 
b-1 



aA;^Ag{t) + Y(;f^)(^ + ^ - «)^^"'^5(«) 



6-1 



E/WaAr^Ag(t) + E 



*^r(M + i) 



(t + M-a)(^-'V(i)5(«)- 



(9) 



Using (6) we get 



6-1 



E/(i)aAr''A5(t) 



6-1 



6-1 t-1 



E /(t)A5W + j;t^ E /W E(^ + /^ - ^(^))^^"'^ Ag(s) 

t—a ^ ^ t—a s—a 

Y^f{t)Ag{t) + —J^J2^9{t) E (5 + A^-'T(t))(^-i)/(s) 



t=a s=a{t) 

6-2 

= fib 1Mb) g{b - 1)] + E Ag{t\A-/^^f{t) 
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where the third equahty follows by Lemma 3.1. We proceed to develop the right 
hand side of the last equality as follows: 

b-2 

f{b - 1Mb) g{b - 1)] + ^ A,9(t),A-^)/(i) 

t=a 

t=b-l ^^^ 

- f{b~lMb)-g{b-l)] + [.g(i),A-^)/(t)]^^^" -^.g'^(i)AGA-^)/(i)) 

t — a 

^/(6-l)g(&)-/(a)g(a)-— -Lg(„) ^ {s + f, - aia)^^-'^ /(s) 

s — (T{a) 
b-2 

where the first equality follows from the usual summation by parts formula. Putting 
this into (9), we get: 

6-1 6-2 

^ f{t)aA^g{t) = f(b l)g{b) f{a)g{a) + ^ (*A^(,)/(t)) g%t) 

t—a t—a 

r(M + i)^ n^^) r(/i + i)^j^^^ 

The theorem is proved. D 

3.2. Necessary Optimality Conditions. We begin to fix two arbitrary real num- 
bers a and /? such that a,/3 S (0, 1]. Further, we put fi = 1 — a and i/ = 1 — /3. 

Let a function L(i, u,v,w) -.T^ x M x K x M — ;► R be given. We assume that the 
second-order partial derivatives L^u, L^v, L^w, Lyy^, L^i, and L^^ exist and are 
continuous. 

Consider the functional £ : J^ — > R defined by 

6-1 

£{y{-)) =Y,my^t),,Atyit)„A^,yit)) (10) 

t—a 

and the problem, that we denote by (P), of minimizing (10) subject to the boundary 
conditions y{a) = A and y{b) — B {A,B e R). Our aim is to derive necessary 
conditions of first and second order for problem (P). 



Definition 3.3. For f E T we define the norm 



11/11= max |/'^(t)|+ max UA^/(i)|+ max |tA^/(t)|. 

A function y E F with y{a) = A and y{b) = B is called a local minimizer for 
problem (P) provided there exists (5 > such that L(jj) < C{y) for all y G 7^ with 
y{a) ~ A and j/(fe) = B and \y — y\\ < S. 

Remark 4. It is easy to see that Definition 3.3 gives a norm in T. Indeed, it 
is clear that ||/|| is nonnegative, and for an arbitrary f £ T and fc G R we have 
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11^/11 = l^ill/ll- The triangle inequality is also easy to prove: 

||/ + .9||=max|/(t)+g(i)|+maxUAn/ + 5)(t)|+max|,A^(/ + 5)(t)| 

< max[|/(i)| + |.g(t)|] +max[UA^/(f)| + |aA?.g(i)|] 

+ max[|,A?/(t)| + |,A?.g(t)|] 

< 11/11 + ll.9i|. 

The only possible doubt is to prove that ||/|| = implies that f{t) = for any 
i S T = {a,a + 1, . . . ,6}. Suppose ||/|| = 0. It follows that 

max|r(f)|=0, (11) 

maxUA^/(f)HO, (12) 

max|,A^/(t)|=0. (13) 

From (11) we conclude that f{t) = for alH S {a + 1, . . . , b}. It remains to prove 
that /(a) = 0. To prove this we use (12) (or (13)). Indeed, from (11) we can write 

aAUit) = A (fxT^i^* + l-a- a(.))(-")/(s) j 

f^(t + 2 - « - a(s))(-")/(5) _ ^(f + 1 - a - ^(5))(-«)/(s) 



^^ {it + 2-a- a(a))(-")/(a) - {t + 1 - a - a(a))(-")/(a)) 



r( 



•^^"^ -A(f + l-a-a(a))(-") 



r(l-a) 

and since by (12) aA"/(i) = 0, one concludes that /(a) = (because {t + 1 — a — 
cr(a))'^^") is not a constant). 

Definition 3.4. A function 77 G J-" is called an admissible variation for problem (P) 
provided r/ ^ and 77(a) = i]{b) = 0. 

The next theorem presents a first order necessary condition for problem (P). 

Theorem 3.5 (The fractional discrete-time Euler-Lagrange equation). If y £ F is 
a local minimizer for problem (P), then 

LM{t) + *A^(,)i„[y](i) + aAf i„[y](t) = (14) 

holds for all i e T'^ , where the operator [■] is defined by 

[y](s) = (s,y'^(s),,A?2/(s),,Afy(s)). 

Proof. Suppose that y{-) is a local minimizer of C{-). Let i]{-) be an arbitrary fixed 

s s 

$(£)=/:(y(.) + £77(-))- (15) 

This function has a minimum at e = 0, so we must have <I>'(0) = 0, i.e., 

b-l 

^ [Lu[y]{tW{t) + LM{t)a^tri{t) + £»[y](t)*Af r?(t)] = 0, 



admissible variation and define the function $ : ( — -n— rw , n fxi i ) — > M by 
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which we may write, equivalently, as 



6-2 



L„ [y] {t)if (t) I t=p(6) +Y.^^ [y] (*)^" (^) 

i— a 

+ 5^L„[y](t)aA,",7(t) + ^L,4y](i)tA^r;(t) =0. (16) 



6-1 



t^a 



Using Theorem 3.2, and the fact that 77(a) = 77(6) = 0, we get for the third term in 
(16) that 



fc-i 



6-2 



Y,LAy]it)aAtrjit) = J2 U^(,)L4y](t) ^f{t) 



(17) 



Using (8) it follows that 



fa-i 



^L^[y](t),Af77(t) 

t—a 

b-1 

= -^L,4y](f)A(,A-''r;(f)) 



i— a 

5-1 



E^-[y]W 



t — a 
/b~l 



A;(^)A77(t) 



r(^ + i) 

vq{b) 



' + u-a{t))^-'-'^riib) 



b-l 



= - ^ L,Jy](i),A;(^) Ar;(t) - — ^ ^(6 + z. - a(t))(''-i)L.4y](0 



(18) 



We now use Lemma 3.1 to get 



b-l 



6-2 



5^L„[y](t),A;(^,)Ar7(i) 

t—a 

6-1 

= ^L.4y](0A7?(t) + — ^^i„,[y](t) ^ (s + ^-a(t))('^-i)A^(s) 

i=a ^'^ ' t=a s=(j(t) 

&-1 6-1 t-1 

6-1 

= 5]A77(t)aAr''i.„[y](i)- 



6-1 



f— a s— a 



(19) 
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We apply again the usual summation by parts formula, this time to (19), to obtain: 

b-l 



6-2 



6-2 



6-2 
1 t=6-l 



= [7y(i)aAr'^L„[y](t)];:^- -^^'^(t)AUAr'^i:„[y](t)) 

t—a 

6-2 



(20) 

Since 77(a) = 77(6) = it follows, from (19) and (20), that 
6-1 6-2 

J2LMit)tA;^,)A^l{t) = -^,7-(i)aAfi„[y](t) 

t—a t^a 

and, after inserting in (18), that 

b-i 6-2 

J2L^mt)tA^bVit) = ^<(t)aAfL„[y](i). (21) 

t—a t^a 

By (17) and (21) we may write (16) as 
6-2 
^ [iMit) + tA"p^,^K[y]{t) + aA^L^mt)] V^t) ^ 0. 

t—a 

Since the values of ri'^{t) are arbitrary for t S T" , the Euler-Lagrange equation 
(14) holds along y. D 

Remark 5. If the initial condition y{a) = ^ is not present (i.e., y{a) is free), we 
can use standard techniques to show that the following supplementary condition 
must be fulfilled: 

/6-1 



^.(a) + w-^ E(* + A'-«)^^~'^^''[y]w 



6-1 



- Y. {t + ^^-c7ia))^^~'^LMit)\+L^ia)^0. (22) 

t=t7{a) J 

Similarly, if y{b) = B is not present (i.e., y{b) is free), the equality 

Luip{b)) + L,{pib)) - L^ip{b)) 

/6-1 
-^ 5](& + 7.-a(t))(''-i)L,„[y](0 



r(i. + i)\^^ 

^ -' \t — a 



6-2 



Y^ipib) + V a{fyt-^^L^\y\(i) = (23) 
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holds. We just note that the first term in (23) arises from the first term on the left 
Irand side of (16). Equalities (22) and (23) are the fractional discrete-time natural 
boundary conditions. 

The next result is a particular case of our Theorem 3.5. 

Corollary 1 (The discrete-time Euler-Lagrange equation - cf., e.g., [12, 14]). If y 
is a solution to the problem 

6-1 

^(y(-)) = E ^(^' y(* + 1)' ^yW) -^ ™n ^^^^ 

y{a)^A, y{b) = B, 
then Luit, y{t + 1), Aj/(t)) - AL„(t, y{t + 1), l^y{t)) ^ for all t e {a, . . . ,b - 2}. 
Proof. Follows from Theorem 3.5 with a = 1 and a L not depending on w. D 

We derive now the second order necessary condition for problem (P), i.e., we 
obtain Legendre's necessary condition for the fractional difference setting. 

Theorem 3.6 (The fractional discrete-time Legendre condition). If y ^ T is a 
local minimizer for problem (P), then the inequality 



Luu[y]{t) + 2L„4y](t) + KM{t) + L,,.v[y]{<t)){^l ~ if 



2 



+ 2(1/ - l)L,^[y]{t) + 2(p - l)L,u.[y]{a{t)) + L,„^[y](t)(l - i/) 

holds for allte T-', where [y]{t) = (t, y-(t), ^Af y(i), tAfy(t)). 

Proof. By the hypothesis of the theorem, and letting $ be as in (15), we get 

$"(0) > (25) 

for an arbitrary admissible variation r](-). Inequality (25) is equivalent to 

6-1 

E [Luu[m){v^{t)f + 2K4mK{t)aA?v{t) + L,Am)iaA?v{t)f 



+2Lu^[y]{tW{t)iKv{t) + 2U4y]{t)aA?7j{t)tA^,7j{t) + L^.^[y]{t)itA^bV{t)? 

2 

Let T gT'^ be arbitrary and define 7^ : T — S' M by 

n(t) = l 1 if * = -(-); 

10 otherwise. 



>0. 



12 
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It follows that 77(a) = rj{b) ~ 0, i.e., 77 is an admissible variation. Using (7) (note 
that 77(a) =0), we get 



6-1 



Y, [L,,u[y]it)iv^it)? + 2L„„[y](i)?7"(i)aA^77(t) + L„Jy](i)(aA^77(i))2 

t — a 



M 



(r) + 2L„ 



6-1 



r(M + i) 
t-i 

') + L,4y](r) 



-5](t + ^-a(.))(^-i)A77(5) 



+2L™[y](i)r?"(i) A7;(i) 

+L„„[y](t) U77(i) + j;^-^ E(* + -" - ^(^))^^"'^ A,7(.) 



E L^Mit) f A77(t) + ^;^-^E(t + A.-a(.))('^-i'A77(.) 



t='^(T-) 



Observe that 



h-l 



E i™[y](t)Ur^E(* + A'-'^W)^'^"''^^(^) +i™('^(T))(-i+AO' 



t=<T(<T(T)) 



6-1 



t=a{T) 

We show next that 
fc-i 



E L™[y](t) A77(t) + — i^E(^+^-'^(^))^''"'^^''(-^) 



E L,v[y]it)(Y7j^^Y.(t + ^^-'^{s)Y''-'^Arjis) 



t=a{a{T)) 



s—a 
6-1 



E i™[y](t)[^^^"'^n:io"^M+* + '^ 



t=a(CT(r)) 

Let t G [cr(cr(r)), 6 - 1] n Z. Then, 

— -^E(*+^-^(^))''~"^'?(^) 



(t - r)r(t - r) 



M 



E(^ + /^-^(s))^''"'^A77(s)+ E {t + ^^-<y{s))^^-'^A^{s) 



r(M + i 

'^ — [(t + M - a(r))(^-i) - (i + A^ - a(<T(r)))(^-i) 



r(M- 



M 



r(M + i) 



r(i + /i-cr(T) + l) 



r(f + ^-cr(cr(r)) + l) 



r(t + M - a(r) + l-i^i-l)) T{t + /i - a(a(r)) + 1 - (^i - 1)) 



FRACTIONAL DIFFERENCE PROBLEMS OF THE CALCULUS OF VARIATIONS 13 



M 






r(f + /i-r) r(t-r + M-i) 



Tit-T + 1) r(t-T) 

(i + yU - T - l)r(f + /i - T - 1) (t - T)r(i - T + ^J - 1) 



(t - T)T{t - t) 

fi in - i)r{t - T + ^i - 1) 



[t - T)T{t - r) 



(26) 



r(/i + i) 


-T — 

-0 


- T)r(i - r) 

^(/x + i + 1) 



(i - r)r(t - r) 

which proves our claim. Observe that we can write tA^77(i) = — t^~(i\^'n{t) since 
77(6) = 0. It is not difficult to see that the following equality holds: 



fo-i 



fo-i 



^2L„„[y](i)7?^(t)tA^r;(i) = - ^ 2L„„[y](i)7?^(i)tA;(^) Ar;(i) = 2U^[y]{T){v-l). 



Moreover, 
&-1 



j22L.4m)aAtv{t)tAUt) 



b-l 



2j2L.^[ym\(Av{t) + Y:r^yY.(* + t^-^{s))^''-'^A,l{s) 



r(M + i) 



fa-i 



Ary(t) + — ^ E (.'^ + ^-^W)^'^-'^A,7(s) 



= 2(i. - l)L„^[y](r) + 2(^ - l)L,^[y]{a{T)). 
Finally, we have that 

b-l 



5]L„.[y](t)GA^y(t))= 



a(r) 



T-1 



6-1 



ArKt) 



r(^ + i) 



- J2 {s + i^~<y{t))^-'^A7is) 



s=a{t) 



^ L.u,u,[y]{t) 



t—a 



6-1 



^^— E (. + ^-a(f))(^-i)Ar;(.) 



r(t. + l) 



s=cr(t) 



+ L„„,[j/](t)(1 - i/)2 + L.u,n,[y]{cr{T)) 



T—1 r 



+ Lu,u,[y]{T){l - lyf + L^w[y]{<y{T)). 
Similarly as we have done in (26), we obtain that 

^^^-^ [(r + ^ - a{t)r-'^ - (a(r) + u ait)^-'^ 
We are done with the proof. 



(a(T)-i)r(a(r)-t) 



D 
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A trivial corollary of our result gives the discrete-time version of Legendre's 
necessary condition. 

Corollary 2 (The discrete-time Legendre condition - cf., e.g., [12, 17]). If y is a 
solution to the problem (24), then 

LuuMt) + 2Luv[y]{t) + LMy]{t) + L,4y]{a{t)) > 

holds for all t G T*^' , where [y\{t) = (t, j/'^(t), Ay(t)). 

Proof. We consider problem (P) with a = 1 and L not depending on w. The choice 
a = 1 implies /z = 0, and the result follows immediately from Theorem 3.6. D 

4. Examples. In this section we present three illustrative examples. The results 
were obtained using the open source Computer Algebra System Maxima.^ All com- 
putations were done running Maxima on an Intel® Core-'"^^2 Duo, CPU of 2.27GHz 
with 3Gb of RAM. Our Maxima definitions are given in Appendix. 

Example 1. Let us consider the following problem: 

b-l 

Ja{y)=Y.{oA?y{t)f—^mm, y{0)=A, y{b) = B . (27) 

t=o 

In this case Theorem 3.6 is trivially satisfied. We obtain the solution y to our 
Euler-Lagrange equation (14) for the case 6 = 2 using the computer algebra system 
Maxima. Using our Maxima package (see the definition of the command extremal 
in Appendix) we do 

Ll:v-2$ 

extremal (LI , 0,2, A, B, alpha, alpha) ; 

to obtain (2 seconds) 

, , 2aB+(a^ -a^ + 2a) A ^ ^ 

Za'' + Z 

For the particular case a = 1 the equality (28) gives y{l) = ^-^, which coincides 
with the solution to the (non-fractional) discrete problem 

1 1 

5^(Ay(i))' = 5](2;(t+l)-2/(t))'^min, y{0)^A, y{2) = B . 
t=o t=o 

Similarly, we can obtain exact formulas of the extremal on bigger intervals (for 
bigger values of b). For example, the solution of problem (27) with 6 = 3 is (35 
seconds) 

^ ^ (6a'^ + 6a) B+(2a^ + 2a^ + 10a^ -2a^ + 12a) A 
y(l) = ^ 



m) = 



3a^ + 6a^ + 15a^ + 12 
12 a^ + 12 a^ -f- 24 a) B + {a^ + a^ + 7 a'^ - a^ + ia"^ + 12 a) A 
6a4-M2a3-K30a2-h24 



http: //maxima. sourceforge.net 
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and the solution of problem (27) with 6 = 4 is (72 seconds) 

3 a'^ + 15 a*' + 57 a^ + 69 a* + 156 a^ -Vlo? + 144 a 



V(X) = 



24 a^ + 72 a^ + 48 a 

I 



e 



B. 



A 



a^ + 5 a"^ + 22 Off' + 32 a^ + 67 a" + 35 a^ + 54 a^ _^ 72 a , 
2/(2) = z A 



i 

24 a^ + 72 a3 + 120 a^ + 72 a 



B. 



y(3) = 



a^ + 6 a* + 30 a^ + 60 a^ + 117 a^ + 150 a^ - \o? ^ 216 o? + 288 a 



c 

72 cc' + 288 OL^ + 792 o? + 576 o? + 864 



a 



c 



B. 



where 



^ = 4 a^ + 24 a'^ + 88 a^ + 120 o? + 196 o? + 144 , 

C = 24 a*^ + 144 a^ + 528 a^ + 720 a^ + 1 176 a^ + 864 . 

Consider now problem (27) with 6 = 4, A = 0, and B ~ \. In Table 1 we show 
the extremal values j/(l), y(2), y(3), and corresponding J^, for some values of a. 
Our numerical results show that the fractional extremal converges to the classical 
(integer order) extremal when a tends to one. This is illustrated in Figure 1. The 
numerical results from Table 1 and Figure 2 show that for this problem the smallest 
value of Jq, a g]0, 1], occur for a = 1 (i.e., the smallest value of Jc^ occurs for the 
classical non- fractional case). 



m 

1.0 - 
0.9 - 
0.8 
0.7 + 
0.6 
0.5 
0.4 - 
0.3 
0.2 -- 
0.1 








+ 
X 



+ 

X 



Figure 1. Extremal y{t) of Example 1 with 5 = 4, A = O, S = 1, 
and different a's (•: a ~ 0.25; x: a = 0.5; +: a = 0.75; *: a = 1). 
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a 


^(1) 


y(2) 


y(3) 


Jc 


0.25 


0.10647146897355 


0.16857982587479 


0.2792657904952 


0.90855653524095 


0.50 


0.20997375328084 


0.35695538057743 


0.54068241469816 


0.67191601049869 


0.75 


0.25543605027861 


0.4702345471038 


0.69508876506414 


0.4246209666969 


1 


0.25 


0.5 


0.75 


0.25 



Table 1. The extremal values y(l), y{2) and y(3) of problem (27) 
with 6 = 4, A = 0, and B = 1 for different a's. 




0.1 0.2 0.3 0.4 0..5 0.6 0.7 0.8 0.9 1.0 



Figure 2. Function Jq of Example 1 with 6 = 4, A = 0, and B = 1. 



Example 2. In this example we generalize problem (27) to 



b-l 



Ja,;9=5]7i(oAry(<) +72(*Af2/(t)^min, y(,0)=A, y(6) = B . (29) 



t=o 



As before, we solve the associated Euler-Lagrange equation (14) for the case b = 2 
with the help of our Maxima package (35 seconds): 

L2: (gaimna[l])*v"2+(gamma[2])*w~2$ 
extremal (L2 , , 2 , A , B , alpha , beta) ; 



m) 



(2 72 ^ + 71 a^ -jia^ + 2jia) A+ (72 /^^ _ ^2 ;9^ + 2 72 /3 + 2 71 a) B 



2 72 /32 + 2 71 a2 + 2 72 + 2 71 

Consider now problem (29) with 71 = 72 = 1, 5 = 2, ^ = 0, S = 1, and /3 = a. 
In Table 2 we show the values of y(l) and Jq := Ja, aivi^)) for some values of a. 
We concluded, numerically, that the fractional extremal y(l) tends to the classical 
(non-fractional) extremal when a tends to one. Differently from Example 1, the 
smallest value of Jq., a e]0, 1], does not occur here for a = 1 (see Figure 3). The 
smallest value of J^, a g]0, 1], occurs for a = 0.61747447161482. 
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a 


y(l) 


Jc 


0.25 


0.22426470588235 


0.96441291360294 


0.50 


0.375 


0.9140625 


0.75 


0.4575 


0.91720703125 


1 


0.5 


1 



Table 2. The extremal y(l) of problem (29) for different values 
of a (71 = 72 = 1, 6 = 2, A = 0, B = 1, and l3 = a). 




0.1 0.2 0.3 0.4 0.5 n.6 0.7 0.8 0.9 1.0 Q' 

Figure 3. Function Ja of Example 2 with 71 = 72 = 1, 6 = 2, 
A = 0, B = 1, and/3 = a. 



Example 3. Our last example is a discrete version of the fractional continuous 
problem [4, Example 2]: 

1 1 
J„ = ^-(oA^y(t))'-y-(i)^min, 2/(0) =0, 2/(2) =0. (30) 

t=o 
The Euler-Lagrange extremal of (30) is easily obtained with our Maxima package 
(4 seconds): 

L3: (l/2)*v-2-u;$ 

extremal (L3 ,0,2,0,0, alpha, beta) ; 

^(1) = ^^- (31) 

a^ + 1 

For the particular case a = 1 the equality (31) gives 2/(1) ~ ■^■, which coincides with 

the solution to the non-fractional discrete problem 

^ i iAyit)f -y^t)^j2l (y(* + 1) - y W)' - y(* + 1) ^ ^^"^ - 

2/(0) = 0, 2/(2) = 0. 

In Table 3 we show the values of 2/(1) and Ja for some a's. As seen in Figure 4, for 
a = 1 one gets the maximum value of Ja, a sjO, 1]. 

5. Conclusion. The discrete-time calculus is a very important tool in practical 
applications and in the modeling of real phenomena. Therefore, it is not a surprise 
that fractional discrete calculus is recently under strong development. Possible 
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a 


y(i) 


J. 


0.25 


0.94117647058824 


-0.47058823529412 


0.50 


0.8 


-0.4 


0.75 


0.64 


-0.32 


1 


0.5 


-0.25 



Table 3. Extremal values j/(l) of (30) for different a's 




Figure 4. Function Ja of Example 3. 

areas of application include the signal processing, where fractional derivatives of a 
discrete-time signal are particularly useful to describe noise processes [25]. 

In this paper we introduce the study of fractional discrete-time problems of the 
calculus of variations or order a, < a < 1, with left and right discrete operators of 
Riemann-Liouville type. For a = 1 we obtain the classical discrete-time results of 
the calculus of variations [19]. Main results of the paper include a fractional sum- 
mation by parts formula (Theorem 3.2), a fractional discrete-time Euler-Lagrange 
equation (Theorem 3.5), transversality conditions (22) and (23), and a fractional 
discrete-time Legendre condition (Theorem 3.6). From the analysis of the results 
obtained from computer experiments, we conclude that when the value of a ap- 
proaches one, the optimal value of the fractional discrete functional converges to 
the optimal value of the classical (non-fractional) discrete problem. On the other 
hand, the value of a for which the functional attains its minimum varies with the 
concrete problem under consideration. 

This research is in its beginning phase, and it will be developed further in the 
future. Indeed, being the first work on fractional difference variational problems, 
much remains to be done. For example, one can extend the present results for 
higher-order problems of the calculus of variations with fractional discrete deriva- 
tives of any order. Moreover, our work also opens new possibilities of research 
for fractional continuous variational problems. In particular, to prove a fractional 
continuous Legendre necessary optimality condition, analogous to the fractional 
discrete result given by Theorem 3.6, is a stimulating open question. One of the 
referees called our attention to the fact that the initial conditions considered in this 
work have ARMA formats. The problem of finding a general formulation leading to 
the specification of the ARMA parameters seems to be also an interesting question. 



Appendix. The following Maxima code implements Theorem 3.5. Examples illus- 
trating the use of our procedure extremal are found in Section 4. 



FRACTIONAL DIFFERENCE PROBLEMS OF THE CALCULUS OF VARIATIONS 19 

kill(all)$ 
ratprint : f alse$ 
simpsum:true$ 
tlimswitch: true$ 

sigma(t) :=t+l$ 

rho(t) :=t-l$ 

rho2(t) :=rho(rho(t))$ 

Delta(exp,t) :=block( def ine(f 12(t) ,exp) , 
return((fl2(sigma(t))-fl2(t))) )$ 

p(x,y) := (gamma (x+1) )/ (gamma (x+l-y) )$ 

SumL(a,t ,nu,exp) :=block( 

define(f l(x) ,exp) , 

f l(t)+nu/gamma(nu+l)*sum((p(t+nu-sigma(r) ,nu-l) ) *f l(r) ,r,(a),((t-l))) 

)$ 

SumR(tl,b,nul,expl) :=block( 

def ine(f2(x) ,expl) , 

f2(tl)+nul/gamma(nul+l)*sum((p(s+nul-sigma(tl) ,nul-l) )*f 2(s) ,s, (tl+1) ,b) 

)$ 

DeltaL (a2 , t2 , alpha2 , exp2) : =block ( 

[alphal :ratsimp(alpha2) ,a:ratsimp(a2) ,t :ratsimp(t2)] , 

def ine(f 3(x) ,exp2) , def ine(q(x) ,SumL(a,x, 1-alphal ,f 3(x) ) ) , 

qO:q(o) , 

o4:f loat(ev (rat simp (Delta (qO,o) ) , nouns) ) , o41 : subst (o=t ,o4) , 

remfunction(f ) , remfunction(q) , return(o41) 

)$ 

DeltaR(t3 , b , alphaS , exp3) : =block ( 

[alphal : rat simp (alphas) ,b:ratsimp(b) ,t :ratsimp(t3)] , 

define(f4(x) ,exp3) , define(ql(o) , (SumR(x,b, 1-alphal ,f4(x) )) ) , 

qlO:ql(z), 

o5:f loat(ev (rat simp (-Delta ( (SumR(x,b, 1-alphal ,f4(x))) ,x)) , nouns) ) , 

o51 : subst (x=t,o5) , remfunction(f ) , remfunction(ql) , return(o51) 

)$ 

EL (exp7 , a, b , alpha7 , beta7) : =block ( 

[a:ratsimp(a) ,b:ratsimp(b) ,alpha:ratsimp(alpha7) ,beta:ratsimp(beta7)] , 

def ine(LL(t,u,v,w) ,exp7) , bl :diff (LL(t ,u,v,w) ,u) , 

sa: subst ( [u=y(sigma(t)) , v=DeltaL(a,t , alpha, y(o) ) , 

w=DeltaR(t,b,beta,y(o))] ,bl) , 
b2:diff (LL(t,u,v,w) ,v) , 
sb: subst ( [t=x,u=y(sigma(x)) , v=DeltaL(a,x,alpha,y (x) ) , 

w=DeltaR(x,b,beta,y(x))] ,b2) , 
sbl :DeltaR(o,rho(b) , alpha, sb) , sbll : subst (o=t ,sbl) , 
b3:diff (LL(t,u,v,w) ,w) , 
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sc: subst ( [t=x,u=y(sigma(x)) , v=DeltaL(a,x,alpha,y (x) ) , 

w=DeltaR(x,b,beta,y(x))] ,b3) , 
sc2 :DeltaL(a,p2,beta,sc) , sc22 : subst (p2=t, sc2) , 
return (sa+sbll+sc22) 
)$ 

ELt(exp8, a, b, alphas, betas, t8) := 

r at s imp ( subst (t=t8, EL (expS, a, b, alphas, bet aS) ) )$ 

extremal (L , a , b , A9 , B9 , alpha9 , beta9) : =block ( 

[a:ratsimp(a) ,b:ratsimp(b) , alpha :ratsimp(alpha9) , 

beta: rat simp (bet a9) ,A1 : rat simp (A9) , Bl : rat simp (B9)] , 

eqs:makelist(ratsimp(ELt(L,a,b,alpha,beta,a+i)) ,i,0,ratsimp((rho2(b)-a))) , 

vars:makelist (y (rat simp (a+i)) , i , 1 , rats imp ( (rho(b)-a) ) ) , 

Xi: [a] ,Xf : [b] ,Yi: [Al] ,Yf : [Bl] , 

X:makelist (rat simp (i) , i, 1, rat simp ( (rho(b)-a) ) ) , X : append (Xi,X,Xf) , 

sols: algsys (subst ( [y(a)=Al ,y(b)=Bl] ,eqs) , vars) , 

Y:makelist(rhs(sols [1] [i] ) ,i, 1 ,ratsimp( (rho(b)-a) ) ) , 

Y:append(Yi,Y,Yf), 

return(makegamma(ratsimp(minfactorial(makef act (sols [1] ) ) ) ) 

)$ 
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